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Phenomenological covariant approach to gravity 
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We covariantly modify the Einstein-Hilbert action such that the modified action perturbatively 
resolves the anomalous rotational velocity curve of the spiral galaxies and gives rise to the Tully- 
Fisher relation, and dynamically generates the cosmological constant. This modification requires 
introducing a single new universal parameter. It requires inclusion of neither dark matter nor dark 
energy. 
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In this paper, we shall consider a Lagrangian density for the space-time geometry in the form of the Ricci scalar 
multiplied to an arbitrary functional of the Riemann tensor and its covariant derivatives. In the first section, we shall 
fix the form of the Lagrangian such that the modified Lagrangian resolves the missing mass problem in the galaxies 
without considering dark matter. This resolution is obtained by requiring that the spherical perturbative solution 
to the modified action for a point- like mass, coincides to the metric suggested by 1-3]. This suggestion takes into 
account the anomalous rotational velocity curves 043 and the Tully-Fisher relation [g. We then shall prove that the 
modified Lagrangian dynamically generates the cosmological constant. All this procedure requires introducing just a 
' ' ■ single new universal parameter. 

m ■ 

<n ; 

I. THE ANOMALOUS ROTATIONAL VELOCITY CURVE 

o 

The space-time geometry around a stationary point-like mass distribution can be presented in the following coor- 
Sh , dinates: 

bJQ 

ds 2 = -A{r)dt 2 + B(r)dr 2 + r 2 dfl 2 . (LI) 

■ 

■ The dynamics of the space-time identifies A(r) and B(r) up to the boundary conditions. In the Einstein-Hilbert 
£N| gravity this metric in the asymptotically flat space-time geometry coincides to the Schwarzschild geometry 

A(r) = -i- = l-^, (1.2) 
B{r) r 

where T) x — 2G e " m , m is the mass, and Gn is the Newton's gravitational constant. Ref. [IHH suggests that in order 
to account for the (almost) flat rotational velocity curve of the spiral galaxies, the space-time geometry around any 

■ point-like distribution of observed matter deviates from the Schwarzschild geometry as follows: 



- r— I 

X 



A(r) = 1 - — + earl lnr + °( e2 ) » ( L3 ) 



r 



5_j 

where e is the systematic parameter of the perturbative series. Such a deviation can be caused by the distribution of 
a kind of matter we could not directly observe, or it can be due to the deviation from the Einstein-Hilbert action. In 
this paper, we adhere to the later possibility. Such a possibility relies within the subject of modified gravities. Many 
studies have been conducted on /(i?)-gravity 9], or /(G)-gravity [To| . However neither of these classes of Lagrangians 
fllj | admits the following perturbative solution: 

ds 2 = -A{r)dt 2 + B(r)dr 2 + r W , (I.4a) 

A(r) = l- — + ea(r) + 0(e 2 ), (I.4b) 
r 

B{r)= ^{l~2eb(r)) (I.4c) 
A(r) 
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where 

i_ 

a(r) = ar^ ln(r) , (1.5) 

where a is a constant parameter. We consider that the following action governs the dynamics of the space-time 
geometry 

S = [ d 4 x^/~detgF[R ljkl ,\7, gilu }R, (1.6) 

where F[Rij k i, V, g^ v \ is a dimension-less scalar constructed out from the Riemann tensor and its covariant derivatives 
and R is the Ricci scalar. This is an extension of f(R) gravity. F must have the following expansion 

F[Rijkl,V, gil A = l + eL[R ljkl ,V,g^}+0(e 2 ) . (1.7) 

For the moment, we do not assume any constraint or form on L[Rij k i, V, g^v}. We consider it as an arbitrary functional 
of the Riemann tensor and its covariant derivatives. We aim to find the simplest form for L[Rij k i, V, g^ v \ for which 
(11.41) is a perturbative solution. In the following, we wish to first demonstrate how we will obtain the perturbative 
equations of motion. 

Let the metric be presented as follows 

g = g (0) + eg (i) + ... : (L8) 
where its indices are not shown but understood. We present the action as follows 

S = J d 4 x(L^[g}+eL^[g} + ---), (1.9) 

where each LW is a functional of metric and its partial derivatives. Inserting the expansion of the metric into the 
expansion for the action we obtain: 



S = I d 4 xL^[g^}+e I d 4 ^ 1 "' 



2 / J4 



d^X 



S 2 g 



•(3 (1) ) 2 



. g M+LV\gW]) 



e=0 



6LW 



Sg 



e=0 



) , (LIO) 



Note that the variation of the first term induces the equation of motion for g(°\ These are the classical equations. We 
fix g(°> such that the classical equations are solved. Having done so the linear term in e in (II.10[) becomes a constant 
number. We then take the variation of the quadratic terms in e to obtain the equation of motion for g( l ) . This process 
gives a non-homogeneous linear equation for g^ 1 '. 

In the following, we apply the above procedure to the following action 

S = J d A x y/- det 5 R (1 + eL[J2y«, Vi,0«]) , (1.11) 



in order to find the equations of motion for perturbation around the Schwarzschild metric (|I.4I) . It is known that the 
functional variation does not generally commute with imposing symmetries on the solution; here we have imposed 
spherical symmetry and time translation. We notice, however, our imposed symmetries are isometries of the Riemann 
manifold -supposedly a smooth manifold- so the principle of symmetric criticality is met [13] • In other words the 
functional variation of the action computed for (|I.4p . indeed gives the equations of motion of A(r) and B(r). 
The non- vanishing components of the Ricci tensor calculated for (II. 4[) read 1 : 

A" A' A' B' A' , . 

R « - 2B-^A + l3 ) + VB> (L12) 
A" A' A' B\ B' , . 

R - = --a + -a^ + ^ + Vb> ^) 



1 Note that Ry, v = R\ x „ 
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where 



1 r ,A! B\ , 

Using the explicit forms of the Ricci tensor, it is straightforward to show that 

-A-v/- det gR = Q + eR {1) + e 2 i? (2) + • • • , (1.15) 
smo 

= -{r 2 a" + 4rb' - 3r h b' + 2a + 46 + W) , (1.16) 

_ -Ar(ab)' -r 2 (a'b)' -4b 2 -2ba+(9r h -l2r)bb' . (1.17) 

One also can use algebraic tools like Grtensorll to calculate the explicit form of RS 1 ' and W- 2 ^ in term of a{r) and 

b{r). Since the Ricci scalar vanishes at the leading order, only the leading term in the e expansion of L[Rijki, V,, gij] 

contributes to the equations of motion of a and b. We therefore replace L[Rijk, Vj, g^] with its value around the 
Schwarzschild metric 

L(r) = L[Ri jk i, Vi,ffy]|e=o • (1.18) 

Doing so the action reads 

S = 4tt / drdi(0 + e£ (1) + e 2 ^^ +•••)(! + eL(r) + 0(e 2 )) (1.19) 



where integration on 9 and <p are performed. The part of the action that contributes to the equations of motion of a 
and b then reads 

S[a,b] = 4ne 2 J dtx J dr(RW + RWl) , (1.20) 
Performing integration by parts, the integration over i?*- 2 -* is simplified to 

/**«-/*<-** + >*> + -., (1.21) 

inserting which into (II.19[) yields 

S[a, b] = 4ire 2 J drdt {-2b' ar + 2b 2 - L(r 2 a" + 4b' r - 3r h b' + 2a + 4b + W)) , (1.22) 

in which the explicit value of M 1 ^ is also inserted . The first variation of 22[) with respect to a(r) and b(r) generates 
the equations of motion for a{r) and b(r); 



<5S[a, 6] 



We first solve (|L24|) for 6(r) 



->■ -2r&' - r 2 L" = , (1.23) 
->■ 2a + 4o + 2ra' + (4r - 3r h )£' = . (1.24) 



b(r)= Cl + ^(L-rL'). (1.25) 



We then insert b(r) into (|I.23|) and solve it for L(r): 

a(r) i 4rci + c 2 

2r 

where ci and C2 are constant of integration. In large r (r > r/,), L can be approximated to 



L = + — — , (1.26) 

1 - ^ ir h - 2r ' V ' 



L « -o(r) - 2 Cl - ^ , (1.27) 
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within which only a(r) depends on the central mass distribution (|I.5[) . c\ and C2 represent some effective cosmological 

parameters in the sense that these constants do not depend on the mass of the galaxy. Around a galaxy, presumably 
we can ignore the cosmological parameters. Therefore we set c\ = = and get 

L(r) « arflnr. (1.28) 

Noting (|TT8|) then yields 

£[%«,Vi )ff y]| e =o « arl lnr > ( L29 ) 

which is an equation that helps us to find the form of Lf-R^jy, V», <7ij]. We would like to find the simplest form for 
£[-Ryfei) Vj, gij] that satisfies (|I.29[) . In so doing we notice that 

GU =Q = (1.30) 

(VG)Vo = 5184 tSr rfe) -518411, (1.31) 

(□ G )U = ^^« 360i, (1.32) 

where G stands for the Gauss-Bonnet term: 

G = R ijk iR ljkl - AR l3 R 13 + R 2 . (1.33) 



So we have 



G \ 12 



(1.34) 



We therefore choose 



(□G)t7 £=0 (360)1 h ' 

h Ug )U=0 w 21nr - (L35) 

G G 

i[%W,Vi,^]| e =o = -a DM a ln(cia* A/ — ) (1.36) 

where a DAf is a constant parameter and ci is a number 2 . We set ci = 1. Also note that a = a DM ( 36 q)3/4 • For this 
choice of L[R ljku ^i,gij] we find 



a 3r^ a a 
2 ~ ~ ~2 



(1.37) 



Note that this value for b means that the rr component of the metric receives no leading correction, this is in agreement 
with the suggestion of [1-3]. To summarize, we have found that 



S= 8^ Jd'xV^TgR (i-^pfrt ln K^)) < 



by construction, admits the following perturbative solution 



(1.38) 



ds 2 = -(1 - ^ + r \ lnr + • • • )dt 2 + (1 + ^ + • • • )dr 2 + rW . (1.39) 

r (360)t h v r 



2 Note that here we need DG in the denominator to account for the Tully-Fisher relation. It is interesting that the existence of derivatives 
in the denominator of the Lagrangian is also required in the phenomenological approach to other problems in Cosmology \\% . 
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Note that the rr component of this metric is not the inverse of the tt component of the metric. 3 A classical (slow mov- 
ing) test probe in this geometry experiences the following effective gravitational acceleration (the leading gravitational 
potential is one half of the tt component of the metric): 

_c 2 r h 12a DM c 2 r h i 

aacc "T^ + ^)F ( ^ )2 ' (L40) 

the first term of which is the Newtonian gravitational acceleration. The second term becomes the dominant term 
when the Newtonian gravitational acceleration is smaller than 

288c 2 2 

Ocritical — (3£Q-)3/2 U dm ■ l 1 - 41 ] 

In the MOND paradigm [TH, [13], the deviation from Newtonian gravitational dynamics is often expected beyond 

mm 

777 

a = (1.2 ±0.27) x 1CT 10 — , (1.42) 



•s 



Inserting this value in the above relation results 

a DM = 1. 78 ^; 2 ? x 10- 13 — L_. (1.43) 

meters 2 

We would like to emphasize that (II.38[) for the above value of a DM can resolve the anomalous rotational curves of the 
galaxies and also reproduce the Tully-Fisher relation [§[. 

We would like to address also the effective gravitational force of (|I.39j) to a relativistic (fast moving) massive object 
(or a probe). The exact effective potential for spacecrafts in the spherical and static geometry of (|I.4a[) is [19] 



1 



-?• 



2 



V ef f(r) = 0, (1.44) 

I 2 E 2 c 2 

Veff ^ r > ~ 2 r 2 B{r) ~ 2 c 2 Mr) B{r) + 2~B{r) ' (L45) 

where dot represents variation with respect to the proper time, and the equatorial plane is chosen to be orthogonal 
to the angular momentum, and E stands for the energy (per unit rest mass) and I represents the magnitude (per unit 
rest mass) of the angular momentum. In an almost flat space-time geometry (A m B « 1) and for a slow moving 
object (E w c 2 ) and in large r, the effective potential can be approximated to 

*®(r)«£ + 3^. C46) 

which is square with saying that the gravitational potential at the leading order is one half of the tt component of the 
metric. When E is large (E > c 2 ) and the space-time geometry is almost flat, the effective potential reads 

Using the expansion series of A and B (|I.4|) . and using (II.37[) . then leads to the following gravitational acceleration 
for a relativistic object 



For non- relativistic object E w c 2 , (|I.48j) converts to (|I.40|) . The dependency on the velocity is due to the fact 
that A(r)B(r) 7^ 1. It implies that the gravitational force of the galaxy (or dark matter halo around the galaxy) 
effectively can be large for massive objects with high speed (v > ). So caution must be taken into account when 
the gravitational acceleration on objects with large speed is being addressed. The stars orbiting the galaxy have low 
speed (^ << 1). So (II.44|) suffices to address the gravitational acceleration they experience. 



3 The corrections to the equations of motion should be such that they do not have vanishing radial null- null component [T^l . 
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II. APPLICATION TO DARK ENERGY 

The cosmological constant in term of h and Q\ reads 

A=Aif2Q A . (II.l) 

Table IV of [20| reports values of J7a and Ho due to +SN la observation as : 

fl A = 0.725±S, (II.2) 

-0.090 
-0.062 ' 

where H a = 100 h km s~ 1 AIPc~ 1 = 3.24 x 10 _18 /is _1 . Using these values, one finds 



h = 0.599™, (II.3) 



1.206+° f 7 i x HT 52 (II.4) 

~°- 073 meters 2 

In the previous section we have shown that a universal parameter, a DM , (|I.43|) resolves the anomalous rotational 
velocity curves of the galaxies. We can utilize the value of a DM to construct a unit for the cosmological constant: 

[AJ = al M = 1.0+° ; s 4 x HT 51 — 1 —. (II.5) 

meters 

which is quite close to the measured value for the cosmological constant. In this section, we wish to modify the action 
for the dynamics of the space-time such that we account both for the anomalous rotational velocity curves of the 
galaxies and dark energy problem. 

Note that we add the cosmological constant to the Einstein- Hilbert action in order to have the following homogeneous 
solution 4 : 

R\ = An, (H.6) 
R = 4A, 

G = \k\ 
□G = , 

where G stands for the Gauss-Bonnet term p. 331) . The modification that we have presented in (|I.38|) diverges for the 
above solution. This divergence is due to the presence of DG in the denominator of the modification. Let us replace 
□G with (□ + kR)G in |Ll8| to obtain: 

S = [ d A x V- dot gR (l-a DM ln(a 4 , — -) ] , (II. 7) 

J V ((D + kR)G)i DM {U + kR)G') y ' 

where k is a non-vanishing number at order one while R is the Ricci scalar. Sufficiently close to a mass distribution 
where we can ignore the effects of the cosmological constant (|II.7j) is the same as (|I.38|) . (|II.7|) . however, admits 
non-trivial homogeneous solutions. 

To find the exact non-trivial homogeneous solutions of (|II.7|) . we must derive its equations of motion for metric. In 
order to take the first variation of the action with respect to the metric (equations of motion), we first rewrite it as 
follows 

S = J d A xyJ- det gC (II.8) 

C = £[i?,G,DG] = R ( l-a DM r ln(a 4 /rn ° — )) . (II.9) 

1 J V ((D + kR)G)i V DM (a + kR)G'J V ' 



The first variation of the action then reads 

'dC r ,dC _ dC 
0R SR+( dG +n 0DG'- 2 



SS = j d 4 x v^del^ ( —SR +(—,+ n—)5G - -C 9ij 5g^ ) , (11.10) 



4 Note that the Cosmological constant is denned by Rij — ^gijR = — A Mcaaurcd gij which leads to (III . 61) where A Moasurcd = A is inferred. 
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The variation of the Ricci scalar reads 

SR = /.',,- V/" + 2E* afl \7 a V l3 5gV , (11.11) 

where E^p is a tensor. The Gauss-Bonnet Lagrangian in four dimensions is topological. So in four dimensions, it 
holds 

j d*x5(Gy/=dtijj) = ^ 5G- ^GgijSgV + 2Ef jafj V tt V = 0, (11.12) 

where E^ a/3 is a tensor. We shall not need the explicit forms of E^ a p and E^ a/3 . We, therefore, don't write them. 
(HT7T2]) leads to 

SG = -/;<,,/</■■' - 2EP ja0 V a V p 5g» . (11.13) 
Inserting (1IL13[) and pi. Ill) in (|II.10|) leads to the equations of motion 

m*« + % ( i + °^g )9 - ~ V 9 « + 2VQV "(i^ - (i + n ^ E ^ = • ( IL14 > 

We would like to find the following solution of the above equations (a dS space-time geometry): 

ds 2 = -(1 - ^r 2 )dt 2 + ^1 - + r 2 (d9 2 + sin 2 Ode?) (11.15) 

1 g-r 

For this solution, the Riemann tensor is covariantly constant: 

VaRijkl = 0. (11.16) 

Therefore any tensor constructed out from the Riemann tensor and its covariant derivatives are covariantly constant. 
Due to this reason, the equations of motion (|II.14|) calculated for (|II.15|) are simplified to 

dC n G8C 1 n , TT1 _ 

dR Rii+ 2dG 9ij ~2 £9ij = °- (IL17) 

Furthermore we note that (|II.15[) holds 
inserting which in (|II.17|) yields 

which can be simplified to 

When the value of A is such that (|II.20|) holds than (|II.15|) is an exact solution of (|II.7j) . After taking the variation 
of the Lagrangian density with respect to R and G and using R = 4A and G — |A 2 , (III.20[) leads to the equation of 
motion for A. In order to succinctly express this equation we express A in term of a new variable defined below 

a A 

A = -™ X A . (11.21) 

4fc 

The equation for A then gives 

v/|fc|cc + 6~3(-l + 51na;) = 0, (11.22) 

This equation has a (real) solution provided that — 1 + 5 In a; <0. So we set the domain of the x variable to a; € [0, ei]. 
In this domain, ([II. 221) can be solved for k; 



Rij = A 5ij (11.18) 
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Also note that we have chosen k to be positive, because x € [0,es] and we want a positive A. Inserting (|II.23|) in 
pi.2ip expresses A solely in term of the x variable: 



V6 



(11.25) 



A = a 4 , v — • (11.24) 

DM 4(1 - 51nx) 2 

We require A = A Moaaurcd while o>dm approximately reads in (|I.43[) . So 

A 1.2 x 10~ 52 

~ (1.78 x 10" 13 ) 4 

DM v ' 

°- 1195 = m V? \s - ( IL26 ) 

4(1 -51ns) 2 

Note that (|II.26|) is solved by x — 0.8892. We, however, want to identify k. The only condition on k is that it should 
be at order one. But we require k to be also either a natural or the inverse of a natural number. According to (|II.23|i . 
x = 0.9167 leads to k = 1. This value of x is not far from the value that solves (|II.26|) . 5 So we set k = 1 and 
x = 0.9167. For these values, the effective cosmological constant reads 

A = 0.1766 a 4 DM . (11.27) 

Now we can read the best value of a DM from (|II.27|) and A = A Mcaaurcd : 

a DM = 1.617t .Sf 2 i x 10~ 13 ^^ , (11.28) 

meters 2 

which is in agreement with (111.51) identified by the critical acceleration in the MOND paradigms. So in summary 

G 



S = [ d 4 x J — det q R \ 1 — a„,, — r ln(a£ 

J \ ((n + R)G)i 



) , (11-29) 



((□ + i?)G)4 v DM (a + R)G 

includes only a single parameter, and dynamically generates the measured value of the cosmological constant, per- 
turbatively resolves the anomalous rotational velocity curves of the galaxies and gives rise to the the Tully-Fisher 
relation. 



III. SUMMARY AND OUTLOOK 



We have shown that a covariant modification to the Einstein-Hilbert action can resolve the missing mass problem 
in galaxies, gives rise to the Tully-Fisher relation, and exactly the same modification dynamically generates the 
cosmological constant. We have constructed the simplest form for the action that is capable of performing these tasks 
pi.29|) . This modification has one single new universal parameter, a DM (|II.28|) . So the suggested modified action 
unifies the missing mass problem in the galaxies and dark energy problem before resolving both of them. 

Our model is not the first geometric model to unify and resolve dark matter and dark energy problems [21H231 ] . 
The previous resolutions, however, do not take into account the Tully-Fisher relation. The metric that they assign to 
a galaxy has a free parameter and that free parameter depends on the baryonic mass of the galaxy. In comparison 
with these, our model dynamically generates the Tully-Fisher relation, and describes the space-time geometry around 
a galaxy with an arbitrary mass only with a universal parameter. Its action, however, seems more complicated than 
the previous models. However, since it needs only one parameter to describe various phenomena we tend to evaluate 
it 'physically simple'. 

The covariant resolution we suggest here, however, should be still counted on as a toy model because the ultimate 
covariant resolution of the dark energy and dark matter problems must address all the phenomena assigned to dark 



matter such as the weak gravitational lensing, the bullet cluster [24|, dynamics of clusters of galaxies, and the CMB 
data. The stability of the cosmological perturbations also must be addressed. Let it be emphasized that, however, 
each of these phenomena demands a set of finite conditions on the form of the Lagrangian governing the dynamics of 
the space-time. Since there exist infinite possibilities available for the form of the Lagrangian, one expects to find a 



5 x = 0.8415 leads to k = 2. x = 0.8892 yields k = 1.3. So k = 1 is the best choice. 
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set of Lagrangians that meets all these conditions. In this paper we already have illustrated that seemingly separated 
phenomena such as cosmological constant and the missing mass problem in galaxies can be due to a single cause. We 
hope that further investigations render that some additional phenomena assigned to dark matter and dark energy can 
also be stemmed from a common single cause: a covariant modified Lagrangian for the dynamics of the space-time. 
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